Theory of high energy features in angle-resolved photo-emission spectra of hole-doped 

cuprates 
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The recent angle-resolved photoemission measurements performed up to binding energies of the 
order of leV reveals a very robust feature: the nodal quasi-particle dispersion breaks up around 
0.3-0.4eV and reappears around 0.6-0.8eV. The intensity map in the energy-momentum space shows 
a waterfall like feature between these two energy scales. We argue and numerically demonstrate 
that these experimental features follow naturally from the strong correlation effects built in the 
familiar t-J model, and reflect the connection between the fermi level and the lower Hubbard band. 
The results were obtained by a mean field theory that effectively projects electrons by quantum 
interference between two bands of fermions instead of binding slave particles. 
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Recently, several groups performed independent mea- 
surements of the electron structure in hole-doped 
cuprates at binding energies E up to leV.[l-5] They 
observed that starting from the nodal Fermi point the 
nodal-direction quasi-particle dispersion breaks up near 
the momentum (tt/A,tt/A) while approaching the zone 
center at E > Ex ~ 0.3 — 0.4 eV. The dispersion curve 
then drops in a waterfall-fashion up to E > E% ~ 0.6—0.8 
eV, where spectral weights reappear while dispersing to- 
ward the zone center. The waterfall also appears in 
the antinodal direction near (tt/2,0). These features are 
observed in the under-, optimal as well as over-doped 
regimes, below or above the superconducting transition 
temperature in hole-doped cuprates. In contrast, this 
feature does not appear in manganites,[4] signifying the 
unique property of cuprates. Given the robust phe- 
nomenology, the mechanism should be independent of 
pairing. Phonons were seen to play important roles at 
low energy scales [6] and at low doping levels [7]. However, 
it is not clear whether they could cause a dynamical gap 
of the order of electron volt. It is also not clear whether 
the polaron physics [8] applies where doped holes are al- 
ready very metallic. The purpose of this paper is to show 
that the robust waterfall feature may reflect the generic 
property of one-band t-J model, serving as a connection 
between the low energy quasi-particles and the residual 
lower Hubbard band (LHB) at higher binding energies 
dominated by local Mottness. In reaching this conclu- 
sion, we deal with the strong correlation effect built in 
the t-J model by projecting electrons with quantum in- 
terference between two bands of fermions. We argue that 
this procedure satisfies the local sum rules for projected 
electrons already at the mean field level, and is there- 
fore able to pick out the higher binding energy degrees 
of freedom. 

We first summarize the main results in comparison 
with the angle-resolved photoemission (ARPES) data. 
Using parameters suitable for hole-doped cuprates, we 
calculated the electronic spectral weight A(h,u>) as a 
function of momentum k and binding energy E = —lu. 
Along the cut (-7T, —it) — ► (0,0) — > (ir,n) — > (ir, 0) — > 
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FIG. 1: (Color online) Intensity maps of A(k,uj) in the 
energy-momentum space (along high symmetry cuts). The 
doping level is (a) x = 10%, (b) x = 20%, (c) x = 30%. The 
intensity scales with the hotness of the color, as shown by the 
color bars. See the text for details. 



(0, 0) this is presented in Figs.l as color intensity plots at 
doping levels x = 10% (a), 20% (b) and 30% (c). At low 
binding energies, the nodal quasi-particle spectral weight 
fades away while approaching the zone center, producing 
a break up roughly at ±(?r/4, 7r/4). On the other hand, in 
the antinodal direction, there is also a break up of the dis- 
persion near (tt/2, 0). The energy scale of the this break 
is around E\ = 0.3el/. The low energy spectral weight 
scales with the doping level, as seen from the increas- 
ing brightness with increasing doping in Figs.l. This is 
the general behavior of doped Mott insulators. At higher 
binding energies E > E 2 ~ 0.6 — 0.8eV, spectral weight 
reappears and behaves in disguise as the missing low en- 
ergy dispersion pushed to higher binding energies. We 
note that in Figs.l the spectral weight in the waterfall 
energy window is weak but nonzero. [9] In order to re- 
veal the waterfall in our case we plot A(k, w) as intensity 
maps for x = 20% in the momentum space at different 
energies in Figs.2(a-e). We also extract a view field from 
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FIG. 2: (Color online) Panels (a)-(e): Intensity maps of 
A(k, ui) in the momentum space for x = 20%. The bind- 
ing energies increases linearly from to 0.8eV. Panel (f): A 
view field extracted from Fig. 1(b) but re-plotted in a log scale 
for the intensity to reveal the waterfall. 



Fig. 1(b) and re-plot it in Fig. 2(f) using a log scale for 
the intensity. The waterfall feature and the pinning of 
the waterfall spectral weight in the momentum space are 
now obvious. In ref.[2] the pinned momentum boundary 
was emphasized as along a linearly-halved magnetic zone 
boundary. However, this does not seem to be a universal 
feature in all cuprates, even though the waterfall indeed 
appears near this boundary. [1, 3-5] The breaking mo- 
mentum is doping dependent as explicitly pointed out in 
ref. [5]. This is also the case in Figs.l. We conclude that 
the qualitative features of our results are in remarkable 
agreement with the data. 

We add that the waterfall energy window decreases 
before it disappears at extremely high doping levels. The 
tendency is seen in Figs.l. In the opposite limit, as x — > 
the high energy features remain (not shown here) and 
this is why we believe that they reflects the LHB. In the 
experimental case, [10] the waterfall seems to also exist in 
this limit, together with the nodal dispersion for binding 
energies larger than but close to 0.356^. The latter may 
be identified as the low energy dispersion in our case 
but folded by the anti-ferromagnetic order, which will be 
studied elsewhere. 

As for the mechanism of the waterfall, in our case it 
is merely a manifestation of a connection between the 
Fermi level and the parent LHB. The higher-energy dis- 
persion is the residual of un-doped systems. By dop- 
ing a Mott insulator, the Fermi level does not have to 
sink immediately within the LHB (the upper Hubbard 
band is pushed to infinity in the t-J model), but remains 
above the parent band. Spectral weights are transferred 
from the parent band to the Fermi level. The amount 
of spectral transfer increases with doping, and is highly 
anisotropic in the momentum space due to the dispersion 



of the parent lower band. In the theory to be described, 
an effective break up in the low binding energy dispersion 
occurs due to the destructive quantum interference effect 
near the zone center, while the spectral weight in the wa- 
terfall energy window arises from the tails of excitations 
near the window boundaries where van Hove singularities 
appear in the Bogoliubov de Gennes bands. One of the 
singularities is already visible around the tips of the high 
energy dispersion in Figs.l. 

The theoretical starting point is the familiar one-band 
t-J model, 
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where no-double-occupancy is implicitly assumed. Here 
n = x, y, x + y, x — y,2x,2y denotes relevant indepen- 
dent bonds of hopping. This model has received tremen- 
dous efforts since the discovery of cuprate superconduc- 
tors. One approach to deal with the projected electrons 
is to regard the electron as a composite particle of a bo- 
son and a fermion.[ll] The slave-boson mean field theory 
(SBMFT) yields a fermion band re-normalized by the 
doping level. The dispersion of quasi-particles in such 
a theory does not break up. One could go beyond the 
mean field level by integrating out internal gauge fields 
coupling to the holons and spinons, but no controllable 
approximation is available for this purpose. On the other 
hand, a similar re-normalized mean field theory (RMFT) 
is obtained by directly resorting to Gutzwiller projecting 
a trial BCS-like wave function. [12] A common feature of 
SBMFT and RMFT is the conventional idea of filling 
the one-band fermion levels. This restricts the extent of 
many-particle entanglement, and predicts that the Fermi 
level sinks immediately within the LHB once a Mott in- 
sulator is doped. We add that there are also efforts to 
bosonize the t-J model, [13] but digging out the fermionic 
excitations turns out to be difficult, except possibly for 
the one-hole problem. [14] 

Recently, a new representation of the t-J model in 
terms of two-band fermions was proposed. [15] The idea 
is to let one of the band, say the p-fermion band, to 
carry the spin and charge of doped particles, while the 
other singly-occupied half-filled band, say the /-fermion 
band, to reflect the neutral spin background. By enforc- 
ing the condition that each p-fermion pairs up with an 
/-fermion into an on-site spin-singlet the allowed Hilbert 
space can be mapped exactly to that of the t-J model. 
In the context of hole-doped cuprates, the p-band carries 
oxygen holes (in the non-bonding band) while the /-band 
carries copper holes, and the p-f spin-singlet is nothing 
but the well-known Zhang-Rice singlet, [16] which maps 
to the vacancy in the effective one-band t-J model. To 
cope with this analogy, we switch from now on to the 
hole picture. 

According to the two-band fermion representation, a 
physical hole removal (or equivalently an electron ere- 
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ation) is such that we annihilate a singlet p-f pair, and 
insert back a /-fermion with the right spin quantum num- 
ber, 

4 = <k = ^| E e aa , ft, {f m f m ), (2) 

where e is a 2 x 2 antisymmetric tensor. No double occu- 
pation of /-fermions is implicitly assumed. We might also 
need to require no double occupancy of the p-fermions. 
But this can be relaxed for two reasons. First, the density 
of p-fermions is low so that the probability of their dou- 
ble occupancy is small. Second, the p-fermions eventually 
tries to form singlet pairs with the /-fermions, which is 
automatically optimized by no double occupancy of p- 
fermions. The hamiltonian can therefore be rewritten as, 

h = E *n(4<w + h - c -) + S I ■ s j e ? e j- ( 3 ) 

{i,n)a (ij) 

Here denotes the spin carried by /-fermions, and 
E p = p-\p\pip\ picks up the site with no p-fcrmions at 
all. Clearly the hamiltonian conserves the /-fermion oc- 
cupancy on each site. [17] As usual, the constrain for /- 
fermions is relaxed to a global lagrangian multiplier at 
the mean field level. The number of p-fcrmions is only 
globally conserved, and can be fixed by a chemical po- 
tential. The total number of p-f spin-singlet pairs is also 
conserved by the hamiltonian, and can therefore be fixed 
by a Lagrangian multiplier. 

We now consider the mean field decoupling of the ki- 
netic term. In the spirit of infinite dimension, in the 
mean field average (rfj cr rfj CT )o we only retain those contri- 
butions with only one inter-site Wick contraction. This is 
equivalent to Wick contract within the <i-operator, leav- 
ing single /- or p-operator. Starting from the identity 

d ° = -j= E e ™' ft' (Apt - ftPi ) 

= E tif" +E/V/-)< ( 4 ) 

* a' a' 

we realize that at the mean field level d a could be ex- 
pressed as a linear superposition, 

d a ~ -L(p ff + 4>U + sY / ^'fl>)- (5) 

Here <f> = f^Pa)o is the mean field p-f particle-hole 
amplitude, and 5 = (pj./f —Pif\)o is the mean field p-f 
pairing amplitude. Indeed 4>f a and 6 e aa i f\, carry 
the same spin and charge quantum numbers as p a does, 
given the fact that f a are charge neutral spinon degrees 
of freedom. 

We observe that the constrain on the number of p-f 
singlet pairs is equivalent to (S* ■ S p ) = — 3x/4, where 
S p is the p-fermion spin, if the average were performed 



exactly by taking into the no-double occupancy of f- 
fermions. In terms of the mean field average we have 
(S* ■ S p ) — 2(Sf ■ S p )o where the re-normalization factor 
of 2 accounts for the effect of projection on /-fermions. 
In a mean field decoupling (with no static spin moments), 
(S f ■ S p )o = -i{4> 2 + S 2 )/8. Wc therefore arrive at an 
important constrain on <f> and 5 as, 

\4\ 2 + \8\ 2 = x. (6) 

This condition must be embedded in the mean field the- 
ory. As one of the important differences to our case, it 
is relaxed in the numerical calculations in ref.[15], where 
5 is set to zero, and (f> — is taken as the signature of 
the spin-charge separated pseudo-gap phase. The latter 
phase is absent due to the constrain in our case. 

From now on wc refer d-operators in the mean field 
sense as in Eq.(5). We argue that combining with the 
above constrain they respect the average but exact local 
sum rules of doped Mott insulators, and as such they work 
as the quasi-particle operator in doped Mott insulators. 
We first recall that in the nonmagnetic uniform states of 
the t-J model, the electron occupied weight is 1 — x, the 
electron unoccupied weight is 2x, and they sum up to 
yield a total weight 1 + a;. [18] These follows from simple 
considerations of the projected electron operators c a c s c\ 
and (? a c s c\. In our case, the total spectral weight of gPs 
is given by 

E(K,4})o = i + H 2 + l<s| 2 = i + z. (7) 

a 

The hole unoccupied weight, or in the reversed picture, 
the electron occupied weight, is given by 

ER4)o = 1 - (x + H 2 + \S\ 2 )/2 = l-x. (8) 

(7 

The electron unoccupied weight is just the difference of 
the above two quantities, 

(dUa)o = (x + 3|0| 2 + 3|«5| 2 )/2 - 2x. (9) 

These sum rules guarantee that we can use the mean field 
Greens function of d-operators even for local operators 
and therefore capture high energy features, with the pro- 
jection arising from quantum interference between two- 
bands of fermions. In contrast, the SBMFT and RMFT 
calculate the electron greens function as G c = G co h+Gi nc 
with an unknown incoherent part. The coherent part 
G co h contains a re-normalization factor given by the dop- 
ing level, and does not satisfy the sum rules alone. The 
condition on (f> and 5 in our case pushes the p-fermion 
to high energies above the Fermi level, and manifests as 
the higher binding energy excitations in the electronic 
ARPES as we demonstrated in Figs.l. 

In the mean field decoupling of the J-term in the hamil- 
tonian, we replace the projection operators E p by (1 —x) 
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for simplest purposes. [15] The remaining decoupling of 
the /-spin exchange is standard, and we arrive at the fol- 
lowing mean field hamiltonian, assuming translation and 
spin-rotational invariance, 

H M f = fk(4 a 3 + A k<?i)fk + Y, p l e k a 3Pk 

k k 

+ ^2\pi^k<J3+V^i)fk+h.c.}, (10) 

k 

where we defined the spinors pk = {pki,P-kl) T 
and/ fe = (/fc T ,/! fei ) T . Here 4 = -(3/4)(l - 

X fj J2n=x,y Xn COS fc„ + (<j) 2 - S 2 ) J2 n *n COS k n - flf IS the 

/-fermion dispersion, Aj. = 2<f)8t n cosfc„ — (3/4)(l — 
x) 2 J X^n=x s cos ^™ * s ^ ne /-fermion pairing function, 
e k = Sn^s.ji cos kn — Hp is the p-fermion dispersion, 
Cfe = Z) n cos k n-^4>, and finally % = £) n <^n cos k n - 
\5. In the above expressions k n — k • r„ with r n the 
vectors along bond n, x« and A„ is the /-fermion hop- 
ping and pairing amplitudes on bond n. Finally fif, fi p 
and A are Lagrangian multipliers that enforce the /- and 
p- fermion occupation and Eq.(6), respectively. Because 
of the constrain on /-fermion occupancy, the effective 
parameters t n and J are re-normalized counterparts of 
their bare values, t n = 2t n and J = 4J, in similar spirit 
to RMFT based on Gutzwiller projection of one-band 
fcrmions.[12] Note that although the renormalized hop- 
ping integral is twice of the bare value the Mott physics 
is built in the quasi-particle-like d-operators as we dis- 
cussed above. Following ref. [15], in the p-fermion disper- 
sion we exclude the nearest neighbor hopping terms in 
the mean field theory. The argument is as follows. The 
bare p-fermions are high energy degrees of freedom, and 
is thus sensitive to local correlations. They would view 
the /-spin background as Neel ordered states, and there- 
fore can only move coherently on the same sublattice. 
This is called coherent path approximation. [15] 

We find that the choice of bare parameters t x , y = 
OAeV, t x+y , x - y = OeV, t 2x ,2 y = OMeV, and J = O.lSeV 
nicely reproduce the experimental features. The general 
conclusion is however not sensitive to parameter tuning. 
Our self-consistent calculations yield that <j> = 5 = \fx~/2 
and [if = 0, and that \ n changes sign while A„ docs 
not in going from rr-bond to y-bond. These combine 
to still form a d-wave superconducting pairing Af? = 
(diidji — d^dji)o- We stress however that the waterfall 
feature is indifferent to this order parameter. For ex- 
ample, the waterfall persists in Fig. 1(c) where A sc — > 0. 
The Matsubara Green's function G^ CT (k, iu n ) for the d- 
operators, which in our case is independent of a, can be 
easily obtained from the mean field hamiltonian. The 
electronic ARPES spectral function is then obtained as 
A(k,u) = -(2/7r)ImGy(k,iw„ -» -w + i0+),[9] where 
the minus sign before u) reflects the fact that we worked 
in the hole-picture in the mean field theory. The re- 



sults are shown in Figs.l and 2, and have been discussed 
previously. More details and results will be presented 
elsewhere. 

To conclude, we proposed a theory that successfully 
explains the waterfall feature in the quasi-particle dis- 
persion observed in hole-doped cuprates. We interpret 
the waterfall as a connection between the fermi level and 
the parent LHB. The theory projects electrons by quan- 
tum interference between two bands of fcrmions. The 
local sum rules satisfied by such a representation make 
it possible to pick up the higher binding energy features 
that is beyond the scope of mean field theories based on 
one-band fcrmions (but of course not beyond the exact 
theory of the one-band t-J model if any). 
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